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We study quasi-normal spectrum of the massive scalar field in the D- dimensional black hole 
background. We found the qualitatively different dependence of the fundamental modes on the field 
mass for D > 6. The behaviour of higher modes is qualitatively the same for all D. Thus for 
some particular values of /xM 1 '' 0- , where fi is the field mass and M is the black hole mass, the 
spectrum has two dominating oscillations with a very long lifetime. We show that the asymptotically 
high overtones do not depend on the field mass. In addition, we present the generalisation of the 
Nollert improvement of the continued fraction technique for the numerical calculation of quasi- 
normal frequencies of D- dimensional black holes. 
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I. INTRODUCTION 

One of the attractive solutions of the hierarchy prob- 
lem is the consideration of our 4-dimensional universe as 
a sub-manifold that is embedded in a higher dimensional 
space-time. The models of type ADD (Arkani-Hamed, 
Dimopoulos and Dvali) |lj allow compactification radius 
of extra dimensions to be of macroscopic size Q- Black 
holes whose horizon radius is much smaller than this 
characteristic length can be well described by solution 
for higher-dimensional static black hole 0. Such black 
holes could be produced at the next-generation particle 
colliders, probably at energies of very low order ~ ITeV 
Q . Thus study of their properties is strongly motivated. 

The very important black hole characteristic is its 
quasi-normal (QN) spectrum: the set of complex fre- 
quencies (see [5| for review). The real parts of the QN 
modes are the frequency of oscillations, while the imag- 
inary parts are their damping rates. These oscillations 
dominate on the intermediately late times of a perturba- 
tion evolution near a black hole. They do not depend on 
the initial perturbation, being defined by the black hole 
parameters only. 

While the QN spectrum of different objects in four- 
dimensional [|| and higher dimensional Q cases was 
extensively studied for massless fields of different spin 
and for the perturbation of the metric itself, our knowl- 
edge about the massive field QN spectrum still has gaps. 
Massive scalar field QN modes were studied only for 
four-dimensional black holes within WKB approximation 
la, and later with continued fraction method of Leaver 
[3 IH1 IHl • Massive vector QNMs were investigated for 
Schwarschild-(anti) de Sitter black holes |12j. Massive 
Dirac QNMs were studied in [l3|. 

It was found that the behaviour of massive fields is 
quite different from that of the massless ones: QN spec- 
trum of massive fields has infinitely long-living modes, 



called quasi-resonances 0, 0] . The late time behaviour 
of the the massive fields in black hole background is the 
same as in pure Minkowski space-time. Asymptotically 
late time behaviour of them does not depend on their 
spin, mass and multipole number [T3 |. Also massive fields 
demonstrate the so-called superradiant instability [l5j . 
which is absent for massless fields. 

This paper is dedicated to the study of the QN spec- 
trum of the massive scalar field in a Z?-dimensional black 
hole background. Note, that the scalar field is widely 
considered as a good model for physical fields when the 
effect of their spin is negligibly small. It appears also 
in cosmological models, and, as one of the fundamental 
states in string theories. 

The scalar field with the mass term can be also inter- 
preted as a self-interacting scalar field within regime of 
small perturbations [l^ . Also, the scalar field gains large 
effective mass when considered in models with extra di- 
mensions of RS-type 01 ■ We study the QN spectrum in 
this regime as well. 

This paper is organized as follows: in Sec.[n]we present 
basic formulae and the method of numerical calculation 
of the quasi-normal modes, Sec. Illll is devoted to our nu- 
merical results and their comparison with those obtained 
in other works. Finally, in Sec. IIVI we summarize the 
obtained results. 



II. BASIC EQUATIONS AND NUMERICAL 
ANALYSIS 



The metric of the Schwarzshild black hole in D- 
dimensions has the form 1181: 



ds 2 = f(r)dt 2 - f^irfdr 2 - r 2 dfl 2 L 



D-2i 



(1) 



where 



D-3 



= 1 - 



16-nGM 
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tj (D-2)Sl D - 2 r D - 3 ' 
Here we used the quantities 
(2tt^ d -^/ 2 

Qd ~ 2 = t((d- i)/2V r (V2) = ^, r(*+i) = *r(z). 
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After separation of angular and time variables, the ra- 
dial part of the massive scalar field equation can be re- 
duced to the form: 



ar ar 



V(r) R(r) = 0, (2) 



where 
V(r) 



1(1 + D -3) 



(D-2)(D-A) 



4? 



f(r) 



D-2 
2r 



f(r) 



I is the overtone number, fj, is the field mass. 

By definition, the quasi-normal modes arc eigenvalues 
of to with the boundary conditions which correspond to 
the outgoing wave at spatial infinity and the ingoing wave 
at the black hole horizon. In 4-dimensional case @, these 
boundary conditions look like 



R 



R 



C + exp(ixrK ro(x+ ^ 2/2x) , as r -> oo , (3) 



C-(r-r y 



'7, 



where x = -\/^ 2 — /J 2 - The sign of \ should be chosen in 
order to remain in the same complex quadrant as to. 

It turns out that the behaviour of the wave function 
in a higher dimensional background differs qualitatively 
from that in the 4-dimensional one. For D > 4 the lead- 
ing term of V(r) — /j 2 at infinity is of order r -2 , that is 
why the appropriate boundary conditions for D > 4 have 
the form: 



R ~ C+ exp(zYj) , as 



oo . 



R ~C_(r-r 



,-iw/(D-3) 



Vh ■ 



(5) 
(6) 



One should note, that these boundary conditions do not 
change the possibility of quasi-resonance existing (see 

M)- 

To solve the equation (J2J , we use the continued fraction 
method |19| | . Equation Q has an irregular singularity at 
spatial infinity r = oo and D — 2 regular singular points: 
r = and the roots of the equation (r/ro) D ~ 3 = 1, one 
of them is the horizon at r = tq . Thus the appropriate 
Frobenius series have the following form: 



R(r) = expfer) (1 



r \ ~ iw /(-D-3) 



where 



i- r -- 



(8) 



is regular in the circle \z\ < 1 for D < 9. For higher 
D, some of the singularities appear in the unit circle and 
one has to continue the Frobenius series through some 



midpoints. This technique was recently developed in |2Cj . 
but in the present work we were limited by D < 9. 

Substituting Q into J2J), one can obtain a (2D — 5)- 
term recurrence relation for the coefficients bi 



min(2D — 6,i) 



E 

3=0 



J2D-5) 



(w) bi-j =0, for i > . (9) 



Each of the functions cj 2£> (u>) has an analytical form 
in terms of s, i and to, for a specific dimensionality D and 
each integer < j < 2D — 6. 

We now decrease the number of terms in the recurrence 
relation 



min(k,i) 
3=0 



(10) 



by one, i. e. we find the Cj / (cj), which satisfy the equation 

min(k— 



J'=0 



(11) 



(4) For i > k, we can rewrite the above expression as 



E c f-M-iM^ = - ( 12 ) 



(fc) / \ 

c fe-M-lM j=l 

Subtracting l|12|) from l|l(J|l we find the relation (|ll|l ex- 
plicitly. Thus we obtain: 

^(uj) = c^ +1 ^(w), for j = 0, or i < k, 

(fc + l)/ N (fc) / \ 

(ft)/ \ 0+1)/ \ C k,i K 10 ) C j-l.i-l{ UJ ) 

This procedure is called Gaussian eliminations, and al- 
lows us to determine the coefficients in the three-term 
recurrence relation 

Co'] &< + cm h-i + c { 2 3 } 6,-2 = 0, for i > 1 
Coi b i + ci 3 i b = 0, 



(13) 



numerically for a given uj up to any finite i. The com- 
plexity of the procedure is linear with respect to i and 
k. 

The requirement that the Frobenius series be conver- 
gent at spatial infinity implies that 



J3W3) ( 3 ) (3) 
n _ (3) L 0.1 L 2.2 Co ; 2 L 2,3 



which can be inverted n times to give 

J3W3) (3 ) (3) (3) (3) 

(3) c 2,n c 0,n-l c 2,Ti-l c Q : n-2 C 2,2 L 0,1 



(14) 



'1,71+1 



,(3) 



„(3) 



c (3) 
c l,l 



r (3) „(3) (3) (3) 
c 0,n+l t -2,n+2 L -0,ra+2 c 2,ra+3 



.0) 

'1,71+2 



.(3) . 
-l,n+3 



(15) 
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The equation l|15|) with infinite continued fraction on the 
right-hand side can be solved numerically by minimising 
the absolute value of the difference between the left- and 
right-hand sides. The equation has an infinite number of 
roots (corresponding to the QN spectrum), but for each 
n, the most stable root is different. In general, we have 
to use the n times inverted equation to find the n-th QN 
mode. The requirement that the continued fraction be 
itself convergent allows us to limit its depth by some large 
value, always ensuring that an increase in this value does 
not change the final results within the desired precision. 

It turns out, that the convergence of the infinite con- 
tinued fraction becomes worse if the imaginary part of 
u increases with respect to the real part. This prob- 
lem was circumvented by H.-P. Nollert [2]]]. To improve 
the convergence of the infinite continued fraction one can 
consider initial approximation for 



b„ 

K-i 



= Rn = 



M 
"2,n+l 



J3) J3) 
L 0,n+l L -2,n+2 



,0) 

'1,71+1 



„(3) 

-1,71+2 



that for large n can be expanded as 

R n {u) = CoH + — M + + 

Jn n 



(16) 



(17) 



(3) 

Since for the four-dimensional case the coefficients c „, 

cj 3 ^ and Cj 3 ^ are known analytically, one can find Cj(u>) 
from the equation 



R n — 



J3) 
-2,ra+l 



.(3) 

'1,71 + 1 



,(3) 



-0,71 + 1^™ + ! 



(18) 



The found expansion l|17|) could be used as an ini- 
tial approximation for the " remaining" infinite continued 
fraction. In order to ensure the convergence of Ijl7|l for 
a given value of u>, one has to start from the found ap- 
proximation deeply enough inside the continued fraction 
ljTo)l. It turns out, that the required depth is less than 
it would be if we start from some arbitrary value as the 
initial approximation. 

As mentioned in |10| . quasi-resonances correspond to 
purely imaginary x- Thus the convergence of the contin- 
ued fraction decreases dramatically, and one needs to use 
the described method for calculation of those frequencies. 

Even though Cg 3 ^, and „ are not known analyti- 
cally in the case under consideration, we are still able to 
find the expansion (|17fl directly from 10. After dividing 
the equation (J5J) by 6i_2r>+6 and using the definition of 
R n = —b n /b n -i, one can find 



2D-6 



2D-7-j 



E (- 1 ) ic S D_8) ( w ) IT ^-fc = ' f°r*»l 



k=0 



(19) 

The important is that the expression (|17fl should be 
defined without any indeterminations. Since the infinite 
continued fraction is convergent, it is clear that R n is 



itself well-defined. Our task is just to recover it from 
the equation 119fl . It turns out that only determination 
of Cq(oj) and Ci(u>) could cause ambiguities: the first 
appears due to a lot of roots of the equation l|19|) and 
the second is connected with two possible signs of y/n in 
(|17|l . Fortunately, both ambiguities can be eliminated by 
using our knowledge about the properties of the series 
||SJ|. Let us consider them more precisely. 

For large n c^^ 5 \uj) oc n 2 , thus Cq(lu) satisfies 



2D-6 



3 J2D-5) 



Hc 



2D-6-J 



{w) = 0. (20) 



3=0 



The equation l|20fl has in general 2D — 6 roots (in fact 
there are repeated roots). One of the roots is always 
Co = — 1 (it is also repeated), implying the unit radius of 
convergence of the series (JSJ . Other roots appear due to 
the existing of additional singular points of the equation 
©. Thus we choose C = -1. 
After fixing Co — — 1 one can find that 



Ci = ±y/-2ir X- 

The sign of C\ can be recovered using the convergence of 
the series (JSJ at spatial infinity (z — 1). Therefore: 

lim b n = 0, 

71 — >00 

i. e. $N : Vn > N \b n \ > |6 n _i|. Since for large n 



b n - 



-Rr, 



-Co 



9l 



we find out that the real part of C\ cannot be negative. 

After the sign of C\ is fixed, the other coefficients in 
1 171) can be found without encountering indeterminations 

Thus we have now possibility to use (|17fl as an initial 
approximation for R n when n 3> roX- This allows us to 
improve the convergence of the infinite continued fraction 
near quasi-resonances. 



III. NUMERICAL RESULTS 

Let us start from an important notion. Since the black 
hole mass in Z?-dimensional space-time cannot be mea- 
sured in units of length, as in four-dimensional space- 
time, we measure all the values in units of the black hole 
horizon radius that is an observable quantity. In order to 
do this we multiply J5J by Tq and consider dimensionless 
quantities [iro and ioro- 

Making use of the technique described in the previous 
section, we have found dependence of the quasi-normal 
modes on the field mass. From the figure^we see that for 
four- and five-dimensional cases increasing of \iro gives 
rise to decreasing of the imaginary part of the funda- 
mental quasi-normal frequency until reaching zero, while 
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Figure 1: Fundamental quasi-normal frequencies (I — 0) for different D and /j,. The frequency for /i = has the largest 
imaginary part. The points were plotted with the step of A/x = 0.1/ro. Solid lines mark the same number of D. 



the real part changes insignificantly. As discussed in 0] j 
when some threshold value of /xro is exceeded, the partic- 
ular quasi-normal mode disappears. The higher thresh- 
old value corresponds to the higher overtone number. 

Since we reached purely real frequency numerically 
with some finite accuracy, we cannot state that purely 
real frequencies exist in the spectrum. We are on posi- 
tion, that at least while /iro tends to some threshold val- 
ues, the quasi-normal ringing lifetime increases infinitely. 

For D > 6 we observe that the imaginary part of the 
fundamental mode tends to zero asymptotically while the 
real part approaches jjL. We observed the same behaviour 
for different multipole numbers. For D = 6 it is shown on 
the figure|21 The picture is similar for higher dimensional 
cases, but since massless field quasi-normal frequency has 
higher imaginary part, one needs to reach higher \xr§ to 
see this tendency. 

Even though the behaviour of the fundamental over- 
tone is different for D > 6, the dependence on /iro of 
higher overtones is qualitatively the same (see figure 
as for D = 4 and D — 5. Also because all the equations 
depend on /z 2 we observe that the contribution of mass for 
higher overtones is of order /i 2 /w as in four-dimensional 
case 0, |nj. Thus the asympotical behaviour of the 
high overtones does not depend on \i. From the figure 01 



one can see that for 7-dimensional case the difference be- 
tween massless and massive /i = r^ 1 cases is insignificant 
already for 20th overtone. According to the analytical 
formula |22| . the spectrum is asymptotically equidistant 
for all .D: 

Ac = l -f(r ) 
and the real part approaches 

Re(uj) = —f'(r ) In. 3 = T Hawking In 3. 

The difference in the mass dependence of the funda- 
mental overtone for D > 6 leads to another remarkable 
fact. Since for high field mass the imaginary part of the 
fundamental overtone tends to zero only asymptotically 
(see figure^, and the imaginary part of higher overtones 
reaches zero for some finite value of \xr§ (see figure 01 , 
there are some values of /itq for which the imaginary 
parts of two overtones are the same. After one of these 
values is reached the overtones can be distinguished only 
by their real part and the quasi-normal ringing has two 
dominant frequencies in its spectrum. 

It turns out that this could happen only for relatively 
large field mass, thus the imaginary part of those fre- 
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Figure 2: Fundamental quasi-normal frequencies for D = 6 as function of /j,ro for I = 0, 1, 2 (represented as dots, rhombuses 
and triangles respectively). The dashed line corresponds to Re(uj) — fi. 



quencies is very small: for example, for the field con- 
figuration of spherical symmetry (I = 0) in D = 6, the 
smallest value of /iro w 5.823, for which there are dom- 
inant modes with u>r w 5.823 and wr n ~ 0.406. Their 
imaginary parts are equivalently small, thus their life- 
time is extremely high. Thus one could observe the su- 
perposition of that two frequencies at late times of the 
quasi-normal ringing. 

One should note that this exotic behaviour was found 
for relatively large field mass for which the linear approx- 
imation (i. e. without taking into consideration the back 
reaction of the field upon the black hole) could not be 
good. 

The fundamental modes of massless scalar field for dif- 
ferent D are presented in table [I] We see that as D 
increases, the real part of the fundamental overtone in- 
creases faster than the imaginary part. Thus fields that 
live in the bulk are better oscillators than brane-localised 
ones [2^, considered in the context of Large Extra Di- 
mension Scenario (see |27| for review). 

We checked the correctness of our results by comparing 
them with those |2^| obtained with the help of the WKB 
method |24| for higher values of the multipole number. 
The WKB method is known to give very accurate results 
in this regime |25j |. 



Table I: Fundamental quasi-normal frequencies of massless 
scalar field (measured in the units of the black hole horizon 
radius) . 



D 


I 


= 


I 


= 1 


I 


= 2 


4 


0.22091 


- 0.20979i 


0.58587 


- 0.19532i 


0.96729 


- 0.19352i 


5 


0.53384 


- 0.38338i 


1.01602 


- 0.36233i 


1.51057 


- 0.35754i 


6 


0.88944 


- 0.53310i 


1.44651 


- 0.50927i 


2.01153 


- 0.50194i 


7 


1.27054 


- 0.66578i 


1.88140 


- 0.64108i 


2.49678 


- 0.63188i 


8 


1.66879 


- 0.78557i 


2.32073 


- 0.76102i 


2.97469 


- 0.75056i 


9 


2.07943 


- 0.89520i 


2.76402 


- 0.87137i 


3.44882 


- 0.86013i 



One should note also, that massless scalar field quasi- 
normal spectrum coincides with that of the metric per- 
turbation of the tensor type for I > 2, which was recently 
studied [2(j. Our results were obtained independently 
and demonstrate complete agreement. 

IV. CONCLUSIONS 

In the present paper we studied quasi-normal spectrum 
of the massive scalar field near the higher dimensional 
black hole within continued fraction method of Leaver 
with improvement of Nollert that was generalised for the 
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Figure 3: Three higher quasi-normal modes (I = 0) for D = 6 and different fj,. The frequency for fi = has the largest imaginary 
part. The points were plotted with the step of Afi = l/(2ro). Solid lines mark the same overtone number. 



case of /c-term recurrence relation for Frobenius series 
coefficients. 

We found that the fundamental modes for D > 6 have 
qualitatively different dependence on the field mass: their 
real part approaches the field mass while the imaginary 
part tends to zero for high values of fj.ro . At the same 
time, the behaviour of higher modes is qualitatively the 
same as for D = 4 and D = 5 cases: the imaginary part 
decreases quickly and the particular mode disappears af- 
ter reaching some threshold value of /itq. 

These cause a very interesting result: for some partic- 
ular mass, the spectrum has two dominating oscillations 
with an identical very long lifetime. 

As in four-dimensional case asymptotically high over- 
tones do not depend on the field mass and satisfy the 



analytical formula for the massless case: 

w« = T Hawking (\n3 - 2nin). 

Also we found that despite the lifetime of the oscilla- 
tions is shorter for larger D, the increasing of the number 
of the extra dimensions causes black holes to be better 
oscillators. 
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